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ABSTRACT 


An ever-present problem in analyzing thermal images for target signatures is the 
influence of atmospheric effects in the signature observed at significant range. The 
compensation for these effects, mostly atmospheric absorption and scattering and path radiance 
requires accurate knowledge of the meteorological parameters for the area involved at the time of 
the measurements. 

Based both on infrared image files taken during the Electro-optical Propagation 
Assessment in Coastal Environments (EOPACE) experiment together with the EOPACE 
environmental data base and on the SeaRad propagation code to generate radiance, a range 
compensation algorithm is proposed in this thesis. Applying SeaRad output adjusted for the sky 
path radiance, an 11 by 11 matrix of the apparent sea temperatures is constructed in which each 
row corresponds to a different zenith angle and therefore range, and each column to a different 
sea apparent black body temperature. By interpolation all sea pixels in the image are range 
compensated. The ship pixels are range compensated by imposing continuity in the sea ship 
interface. The magnitude of scene temperature correction required is of the order of -2.3+ 1.7 °C 


which is comparable to the precision of the recorded data. 
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P INTRODUCTION 


All predictive codes are based on the target-plane target signature. This poses a 
challenge to range-compensate target signatures for the influence of atmospheric effects and other 
sources of radiation. 

The Electro-optical Propagation Assessment in Coastal Environments (EOPACE) 
database allows an evaluation of this compensation method. A very comprehensive 
characterization of the meteorology and sea surface parameters, including data from two Naval 
Postgraduate School (NPS) buoys moored in San Diego Bay, 1s available for the period of 
intensive measurements in EOPACE. Full shipboard location and local meteorology data were 
recorded on the R/V POINT SUR for two days of operations in the sector south of Point Loma, 
where NPS polarization imaging equipment was located. During these two days the POINT SUR 
operated in close liaison with the NPS-NACIT imaging group on the south tip of Point Loma. 
The POINT SUR took up stationary positions at varied aspect angle at selected ranges along 
several bearing lines from the images. 

The purpose of this thesis was to develop a method to compensate the imagery for the 
atmospheric extinction and radiance effects. But, to check the correctness of the method more 
than just the thermal image data 1s necessary. The NPS Boundary Layer group played that role, 
supplying meteorological and geographical data. 

A usual approach to range-compensate a thermal image is to begin by finding the 
equivalent radiance, here called the “full range radiance”, for each pixel value in the image plane 
(that is, for each apparent temperature). With this full range radiance and calculating the 
atmospheric transmittance and path radiances using SeaRad (a computer code that calculates 
transmittance and path and sea radiance), the radiance at the source plane (zero range radiance) 
can be calculated and also the corresponding zero range apparent temperature (the range- 
compensated temperature). SeaRad is a rigorous geometric capillary wave model based on the 
Cox and Munk wave-slope statistical model suitable for moderate wind, high altitude, slant-path 
sensor configurations, that was developed by Naval Command Control and Ocean Surveillance 
Center (NCCOSC), now Spawar Systems Center, San Diego. 

But this approach is very arduous and computer time-consuming and in consequence a 
new method is presented in this thesis. Basically, using SeaRad, a matrix of apparent sea 
temperatures at full range with dimensions eleven by eleven is built first, in which each column 


corresponds to different arbitrary apparent sea temperatures at zero range and each row 


corresponds to different arbitrary zenith angle. Then, for each row of sea pixels at the image 
plane the corresponding zenith angle is calculated. Using two rows of that matrix, a new matrix 
with dimensions one by eleven is interpolated for that zenith angle, but which keeps the same 
column correspondence of the previous matrix. Finally, for each sea pixel value (that is, each 
apparent temperature at full range), the apparent temperature at zero range can be found by 
interpolation, using the one by eleven matrix. For the ship the process is the same except that the 
base matrix has dimensions one by eleven because all ship pixels have the same zenith angle. 

This thesis begins by first finding the calibration curve for the AGA equipment and for 
the conditions of the time of the experiment. With this curve, all image data were mapped from a 
digital level unit, that is characteristic of the equipment, to temperature. This image represents 
an apparent temperature resulting from the interaction of different sources of radiation. | 

The interactions of radiation from the various sources are explored in Chapter III. 

In Chapter IV, the available input data are presented and analyzed. Also, other necessary 
inputs are determined. Those are the horizon position, the sun position, the target position, etc. 

The core of this work is in Chapter V. First, the structure and the mathematical treatment 
applied to the data in this thesis are presented. Second, the SeaRad results are adjusted for the 
sky path radiance by making the calculated SeaRad apparent temperature at the horizon zenith 
angle to be equal to the measured horizon apparent temperature. It is found that no SeaRad 
correction is necessary for the sea other than the sky path correction. Finally, the range 
compensation method is introduced. 


The conclusions and recommendations are presented in Chapter VI. 


Il. EOQPACE DATA 


In the Electro-optical Propagation Assessment in Coastal Environments (EOPACE) 
experiment conducted on April 9 and April 10 of 1996 in San Diego, a set of polarized and 
unpolarized thermal images were taken. The Naval Postgraduate School (NPS) Boundary Layer 
Meteorology Group and the NPS Center for Infrared Technology (NACIT) also have contributed 
to that specific EOPACE experiment by recording meteorological and geographical data. In 
[Ref.7] and [Ref. 8] there is a complete description of the experiment with characterization of the 
experiment location, recording format and so on. 

The experiment can be represented diagrammatically by Figure 2. 1, where an observer 
(sensor) receives an image attenuated during its propagation through the atmosphere. The 
observer has also available recorded meteorological and geographical data that would allow him 


to analyze thermal vision images. 





Meteorological and 
Geographical Data 





3 Thermal Image 





Figure 2. 1 The EOPACE Experiment. 


In the experiment of April 9 and April 10 of 1996 the following parameters were kept 


invariant: 


2d 


: The observer geographic coordinates (also called sensor or camera position in 


this document): 32° 39’ 55” N, 117° 14’ 31”. 


3 The instrumentation: AGA 780 Thermovision System. 

e The target: R/V POINT SUR: a 135 foot long ship. 

The variant parameters were: 

© The ship position relative to the observer. 

e The ship heading. 

e The uncontrollable ship signature and sea temperature. 

e The Sun position. 

e The uncontrollable weather conditions: wind intensity and direction, relative air 


humidity, air temperature, and air pressure. 
A. INVARIANT PARAMETERS 


The fact that the same AGA 780 Thermovision System equipment was used during the 
whole experiment does not mean that the instrumentation is not adding any sort of additional 
error to the measurements. Two elements can have a great influence on the results: the 
equipment calibration and the equipment drift. Another kind of error that could be significant is 


operator error. 


1. Calibration 

The use of different controller units does not affect the results of thermal image 
acquisition, but different scanners have different calibration curves and will affect the output. In 
the experiment the scanner used was: AGA serial number 5016, locally called “the NOSC head”. 

Since there was no calibration curve directly associated with the data taken in this 
experiment, it is necessary to determine experimentally a calibration curve and check the 
differences in relation to a past determination. Then, the calibration curve found should be 
applied to the experimental data in order to figure out the magnitude of error, if it exists. 

The comparison calibration curve used was one made by E.C. Crittenden on March 21, 
1996 on the same scanner. 

For the calibration determination of this thesis the following settings and equipment were 
used: 

e Scanner serial number 5016. 


° Acquisition software: PTRWIN. 


” Black body standard source Graseby Infrared, Inc. model 140/201 

‘ 7° FOV lens; 

e The Black Body distance from the lens is | meter. 

The PTRWIN program [Ref. 12] is an acquisition and analysis software package that 
uses a 12 bit Analog-to-digital converter for converting the detector signal intensity into a real- 
time image on a computer screen and for digital image recording. PTRWIN has two functions 
that are very useful in analyzing images: the spot and the area options. The spot option of the 
PTRWIN software allows measuring spots of a displayed image, and the area option to calculate 
the statistics for selected areas of the screen. 

The following procedure was adopted for reading the AGA temperature corresponding to 
the Black body temperature: 

@ Use the spot and area options of the PIRWIN software to read the correct value 

of the output in isothermal units (IU). 

° Adjust the thermal level dial until saturation of the highest temperature color 
occurs. At this time the thermal level value presented at the spot camera dialog 
window will fe the same as the maximum IU value displayed in the color bar 
representation of the image (saturation). From that point on, it will stay constant 
if the thermal level is decreased, but the IU value will also decrease, showing that 


the recorded value was the maximum. 


a. Black Body Characterization 
The black body used in the calibration experiment conducted in the laboratory 


has the following specifications: 


® Temperature range: ambient to 230 °C. 

® Absolute calibration accuracy: + 1.5 °C. 

e Stability: + 1 °C. 

> Resolution: | °C. 

2 Cavity opening: 12”x12” square. 

2 Temperature distribution: approximately 2 °C edge to edge. 
e Cavity emissivity: 0.95 + 2%. 


A flat-top profile of temperature along the vertical and horizontal axes of the 


black body cavity is desired. To verify this surface uniformity, the temperature distribution along 


the vertical axis (1) and along the horizontal axis (2) of the black body at 80°C are shown in 


fneune 2. 2. 
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Figure 2. 2 Temperature distribution on the Black Body Cavity at 80 °C. 


Note that the horizontal axis has better uniformity than the vertical. Also note 
that the maximum and minimum values for both curves are not far apart, confirming the 


uniformity of the black body. 


b. Calibration Determination 
Data for three different filters were measured. Filter 6 is unpolarized, filter 5 is 


vertically polarized and filter 4 is horizontally polarized. Table 2. 1 presents those data for an 


ambient temperature of 17 °C. 


found: 


Black Body Filter 
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Table 2. 1 Calibration Data for Ambient Temperature of 17 °C. 


Usually the calibration curve is described by the curve-fit equation: 


A 


i aa ae (rls) 
Ce cia 
Where: 
‘ T 1s degrees Kelvin. 
° A, B and C are constants that are found from curve fitting. 


Using Mathcad 8 software, the following constants for the unpolarized case were 


e A = 418.751 
e B=1115.0 
° C=0.216 


And, for the polarized cases (filter 5 and 4): 


© A = 5.679 
. B =57.46 
e C=0.918 


To see the influence of ambient temperature on the isothermal values for the 
unpolarized filter, another set of measurements was taken at an ambient temperature of 23.4 °C 


with the results presented in Table 2. 2. No significant discrepancies are noted. 
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Table 2.2 Calibration Data for Ambient Temperature of 23.4 °C. 


The calibration constants for the average IU values are: 


@ A=2318 
® B=1218 
a C = 0.742 
ron Previous NPS Calibration 


Previous NPS calibration of the AGA-780 over the time interval from 1990 to 
1996 were made mostly by Dr. E. C. Crittenden and were expressed in terms of the curve fit 
equation incorporated in the CATS acquisition and processing software (see Equation 2.2). 


These calibrations are referred to in this thesis as the Crittenden calibrations. 
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Where: OS is offset. 


The calibration constants for Equation 2.2 and the corresponding values for 


Equation 2.1 are reported in Table 2. 3: 
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Table 2. 3 Crittenden Calibration Constants. 


















d. Comparison Between Crittenden Calibration and New Calibration 

For comparison purposes, in Figure 2. 4 and Figure 2. 5, the Crittenden 
calibration was plotted together with the new calibration for the polarized and unpolarized cases 
respectively. The data agree in shape but there is a non-constant offset between the data. This 
offset is a characteristic of the scanner controller and can be adjusted by rotating the thermal level 


adjustment screw located on the front chassis (number 11 in Figure 2. 3). 
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Figure 2.3 AGA-780 Controller Front Panel. 
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Figure 2. 4 Calibration Comparison for Polarized Filters. 
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Figure 2. 5 Calibration Comparison for Unpolarized Filter. 


This offset would not be serious for analysis of field data if there were 
calibration- checking data (a black body at a known temperature) in the recorded scene data that 
could allow determining the correct calibration curve. This calibration determination could be 
done, first, by finding the isothermal unit IU,.- corresponding to the calibration-check black body 
temperature T,e-. Second, in the laboratory try to minimize the offset for that particular point 
(User , Tre) by rotating the screw located on the front panel of the controller (number 11 in Figure 
2.3). Finally, find the calibration parameters as explained in this subsection. 

Note that to assume that one curve is just the offset of the other can result in high 
error. For example, assume that an offset between the two curves were found for 19.19°C. Also 
assume that the correct curve 1s the new calibration. The offsetting of the entire Crittenden 
calibration curve by a constant amount in order to match its curve with the new calibration at 
19.19 °C is represented in Table 2. 4 and Table 2.5. In the same tables the difference in 
temperature for each IU is presented. These errors in the temperature scale can be seen in Table 


2.4 and Table 2. 5. 
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Figure 2. 6 Temperature Error Using Offset Method. 
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Table 2. 4 Temperature Error for Polarized Filters - Offset Method. 
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Table 2. 5 Temperature Error for Unpolarized Filter - Offset Method. 










Then, by first finding a better calibration curve as described in the beginning of 
this sub-subsection and then offsetting it (if another adjustment is needed), the total error can be 
minimized. 

But a quick look over the files of the experiment will show that there is only one 
black body measurement recorded, which used filter 4 and was made at the end of the day 10 
April 1996. In consequence, although still not ideal, this thesis will be based only on this April 
10, 1996 data, assuming that neither an operational error nor an equipment drift happened in the 


proximity of the time of the black body measurement. 


é. Data for Calibration Checking 

The data measurement for calibration checking was recorded on the file 
“NSFA1093.ptw”, using horizontal polarization, and the black body set to 40°C. A check of the 
data quality is needed knowing that the better the “flatness” of the black body the better the 
quality of the measurement. The method depends on the width of the feature being greater than 
that required for the response to rise to its maximum value. In Figure 2. 7 the temperature 
profiles for the horizontal plane (1) and for the vertical plane (2) across the black body are shown. 
The black body response is not as flat in the horizontal plane showing that a statistical evaluation 
over a Specific area of the black body surface is needed. This can be done using the area option 


of PTRWIN. 
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gure 2. 7 Quality of Black Body Measurement. 
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Applying the area function of the PTRWIN software as shown in Figure 2. 8, the 
IU value corresponding to a temperature of 40 °C is 56.524 + 0.148 IU. 
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Figure 2. 8 Isothermal Unit Value for a Black Body at 40 °C. 


fi Digital Level to Isothermal Unit 
All EOPACE images were recorded in the digital level format. In consequence, 


for further image analysis it is necessary to determine the relation between the digital level and 


the isothermal unit value. 
The digital level (D) ranges from 0 to 4095 (12 bits scale) and the isothermal unit 


from | 7. = to | TL = , where TL is the thermal level and TR is the thermal range. 


Both data are recorded in the image headers. 
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Applying linear correspondence: 


TU 2p ee Ce 
2 4,095 
g. The Correct Determination of Calibration Curve 


Knowing (User , Tres) = (56.524 + 0.148 IU, 40 °C ), and setting a black body at 
40 °C, it should be easy to adjust the controller to provide an output of = 56.524 by rotating the 
screw located on the front panel of the controller, and then find the correct calibration curve. 
Fortunately, the (56.524 + 0.148 IU, 40 °C ) point matches with the data collected on Table 2. 1, 
and in consequence Table 2. | contains the points of the correct calibration curve. Also, 
apparently no big changes have occurred in the AGA system since the end of April 1996 
EOPACE experiment. | 

The relationship of the IU values to apparent temperature can be expressed in the 


calibration equation used in the PTRWIN software, namely: 


EeR alae 


lV = z Z (2. 4) 
eT-F eV —F 
Where: 
° T is the object temperature in degrees Kelvin. 
3 € is the object emissivity. 
e V is the background temperature in degrees Kelvin. 
e R, B and F are constants that are found from curve fitting. 


The constants in the PTIRWIN equation can be related directly to the CATS 
calibration constants as follows. 
In the case in which the scanner field of view is smaller than the black body 


surface, no background is seen by the scanner and V = T in Equation 2.4, thus: 


R 


— Rk F 
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Comparing Equation 2.1 and Equation 2.5 for the polarized case: 


i a = 1.0893 (2. 6) 
C 
R=AeF=6.186 Qa 
And for the unpolarized case: 
Fe eS = 4.630 (2. 8) 
C 
R=AeF =1,938.662 (2. 9) 


Using the above calibration constants in PTRWIN and applying the results to the 
file “NSFA1093.ptw”, a temperature of 42.8 + 0.4 °C can be found for the apparent black body 
temperature. Again, the area function of the PIRWIN software is used as shown in Figure 2. 9. 
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Figure 2. 9 Black Body Temperature Check. 


It may be inferred that the error associated with the polarized calibration curve is 
+2.8+ 0.4 °C. 
Since there is no black body measurement for the unpolarized case on the 


experiment date, some black body measurements recorded in the preparation phase of the 
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experiment were used. All these data were recorded on April 6, 1996, using a black body at 40 


°C and filter 6. The results are tabulated in Table 2. 6. 


File Name Measured Black Body 
Temperature in °C 


NSCA0307.ptw 35.1402 


NSCA0314.ptw 38.4401 
NSCA0339.ptw 39.2402 


NSCA0342.ptw 38.1401 


Table 2. 6 Black Body Temperature for Filter 6. 





Then, the average error associated with the measured unpolarized calibration 


curve can be found using Table 2. 6 to be -2.3 + 1.7 °C. 


h. Influence of Temperature 

In the above calibration curves, the temperature inside the scanner was not 
recorded, but it will expected to be near the ambient temperature, since the experiment can be 
considered as a case of steady state. But to see the effect of temperature on the calibration 
parameters extreme conditions are needed. The conditions required were determined 
experimentally in the drift test described in the Subsection 2 that follows. 

Note that the calibration curves measured in the laboratory for the three widely 
different temperatures are shown plotted together in Figure 2. 10. Although these wide variations 
in the laboratory produce observable changes in the calibration curve, such temperature variations 
would not occur in field measurements. So the calibration variation (Figure 2.5) cannot be 


attributed to the ambient temperature. 
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Figure 2. 10 The Temperature Effect on Calibration Curves. 


a Drift 
To check the drift of the system, the output of a scanner radiated with a 195 °C source 


was monitored as depicted in Figure 2. 11. 
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Figure 2. 11 Layout of the Drift Experiment. 
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Because the black body model 140/201 can be set to a temperature that ranges from zero 
to 230 °C, one black body was used just as a heater source. 

During the whole experiment the target temperature was kept constant at 20°C, and a 
thermistor was installed inside the chassis of the dual scanner, between the two dewars. 

The results are contained on Table 2. 7. The drift test was stopped when the temperature 
inside the scanner dropped with the addition of Liquid Nitrogen (LN), because at this condition 


the system was in steady state. 
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Table 2. 7 Drift Test Results: Observed IU Values Versus System Internal Temperature. 
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Note that the drift 1s significant when the LN level in the dewer is low. This 1s an 
operator error that could easily happen in a field experiment. 

Also the maximum drifting error for this extreme ambient condition was only 2 IU. In 
terms of temperature this corresponds to approximately 2°C, considering that the variation of IU 
ACU) _ 
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Using this experimental apparatus, a calibration data set was collected for this extreme 


is related to the variation of temperature by: 


ambient condition, and 1s shown in Table 2. 8. This calibration curve can be used to check the 


influence of temperature on the calibration constants (see Figure 2. 10). 
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Table 2. 8 Calibration Data for a Temperature Inside the Scanner of ~ 30 °C. 
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3} Operator Error 

All the images were recorded with the AGA scanner equipped with the 3.5° lens. This 
field of view of 3.5° is input manually by the operator and is recorded in the file header of the 
images. Even though it does not modify the measurement results it does influence the choice of 


the correct calibration curve and it effects the geometry of the image seen by the observer. 


a. Influence of Lens FOV on Calibration 
To see the magnitude of the error of choosing the wrong calibration curve, the 
calibration performed on contract by Al Fudali in April 1990 for the scanner Serial Number 401 1 


and for two different lenses have been plotted on the same scale in Figure 2. 12. 
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Figure 2. 12 Calibration Curves for Two Different Lens. 
The error in assuming one field of view in place of another is less than 0.5 IU in 
the region of temperature below 60 °C, as seen in Figure 2. 13. In consequence the calibration 


curve for a 7° lens can be used instead of 3.5° lens without perceptible error. However the pixel 


dimensions will be considerably different. 
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Figure 2. 13 Difference in the Calibration Values for Two Different Lenses 
(2-5 um, and 8-12 tum). 
b. Influence of Lens FOV on Geometry 
The main effect of the field-of-view (FOV) on geometry is in the image size and 
the relation of the pixel locations to range (FOV is the angle in Fig 2.15). Because the range will 
impact greatly the attenuation of radiation throughout the atmosphere, it is necessary to provide a 
procedure to verify the correct FOV. Since what is needed is just an estimate, the earth can be 


assumed spherical and we can adopt the geometry of Figure 2. 14. 
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Figure 2. 14 Range Corresponding to the Down Look Angle to Sea Surface from Sensor 
Elevation. 
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From the point of view of the sensor, it sees the geometry of Figure 2. 15. 
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Figure 2. 15 Image Plane Geometry. 


From Figure 2. 15: 





a ian-( = e fan For (2. 10) 
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And from Figure 2. 14 we obtain the down look angle: 


= re 
E = COS ( (2211) 
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The point “T” , coordinates (xr, yr ) in Figure 2.14, is the solution of the 


following system: 


a. y? = re* 
(Quay 


y=(reth)—xetan(e+<') 


Thus, using Mathcad symbolic menu: 


eae (n+ re)e tan(e + €')— re? eran(e +e’) ~2ereeh—h’ (2. 13) 
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For small h, x7 is a good estimate for the distance from the observer to the 
target, and Equation 2.13 can be used in comparison with the geographical data collected and to . 


check the correctness of the FOV information. 
B. VARIANT DATA 


Meteorological and geographical data were recorded by the NPS Boundary Layer 
Meteorology Group and the NPS Center for Infrared Technology (NACIT). 

Sixteen thermistors were installed on the target ship at the locations shown in 
Figure 2. 16. The file “skindata.pte” contains those data for April 9 and 10, 1996. The data were 
sequentially recorded at a 20 seconds rate, where the first datum is the time in seconds in Pacific 
Standard Time, followed by the 16 measured temperatures from the 16 thermistors. There are 
two blocks of data with time varying from 0 to 86,400 seconds. The first block corresponds to 
April 9 and the second block to April 10. 
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Figure 2. 16 Location of Thermistors. 
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Reference 25.0 C 
Port Below Stack 
Port Bow 
Starboard Stack 
Aur 

Port Stack 

Bow Boathouse 
Starboard Boathouse 


The remaining data are recorded on the file “eop30sec.101”, usually referred to as a 


METOC data file. This file holds only data from April 10 and contains 25 columns of distinct 


information. The METOC data needed for this thesis are those included in Table 2. 9. 


Z 





04 
05 
10 





i 
12 
13 
14 
15 
21 

2 









Column 


Number 


a 
i aa 
a 
a 
_ a 
Se a 
a 


a Ship Latitude, degrees-minutes Chet is 
Ship Longitude, degrees-minutes eam 9 








Description 


Sr a 










Sea Surface Temperature, °C tt.t 


Table 2. 9 Experiment Meteorological and Geographical Data. 
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Wl. INFRARED RADIATION 


From direct viewing of a thermal image, it is hard to identify the modification that the 
radiance of each different source had suffered on its path to the sensor (camera). Taking the 
Electro-optical Propagation Assessment in Coastal Environments (EOPACE) images as an 
example, they have as sources of radiation: the sun, the sky vault, the sea, the target (ship), and 
the atmosphere. The radiation leaving each of these sources can be attenuated, reflected, or 
scattered by other sources. Also the attenuation depends on the distance that the radiation had 
traversed and is also a function of the atmosphere model. The sun component depends on the sun 
position relative to the sensor, the sun reflected on the sea depends on the sea wavelets, and other 


complexities add up to make the range compensation of a thermal image very arduous. 
A. TERMINOLOGY 


The radiation leaving a small area element of source can be represented as in Figure 3. 1 


SQ Op 


ee. 


TRE. 


. 
\ 

”, 

in » 

he e 
a + 
« 
7 


dAe cosO 
Figure 3. 1 Radiation Leaving an Area Element. 


The following terminology can be defined [Ref. 1] and [Ref. 3]: 


e Radiant Flux = 0¢ 9 — watts 


OP y watts 
2 
0A ecos@ ATK 





e Radiant Exitance or Irradiance = Mg = 


6¢, — watts 
oe —» ee 


e Radiant Intensity = Jpg = a 
sr 


a 


e Radiance = Lg -_%9 ____ (Sto, _ watts 
dAecosO' d6AecosOedQ spec? 
Usually the above relations are dependent on the wavelength, and a 1/um factor could 
appear multiplying those units, which then represent “spectral” quantities. 
If A is part of a “point source”, then the radiant flux will be independent of the direction. 
If dA is part of an extended Lambertian source, then the radiant flux sill not be 


distributed uniformly with direction as shown in Figure 3. 2, and the following relation will be 


true: 0¢g = 09, ¢ cos 0 





Oo fe Normal Flux 


Figure 3. 2 Radiant Flux in a Lambertian Extended Source. 
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Then, for a Lambertian source: Lg = ————++—__ = —-—_ = constant, and the 


bAecos9edQ 6A0eSQ 


radiance is independent of the observing direction. 
B. RADIATION INTERACTION 


For targets located in the sea, the radiation that goes into the sensor can be represented as 
depicted in Figure 3. 3, where the specularly reflected sun component refers to the sun radiation 
reflected from the sea and “t” refers to transmittance. The use of Figure 3. 3, Figure 3. 4, and 
Figure 3. 5 will allow us to write the total energy that goes into the sensor and therefore to predict 


the temperature. 
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Figure 3. 3 Sources of Radiation for the EOPACE Data. 


The sun depicted in Figure 3. 3 1s not in the field of view (FOV) of the sensor. However, 
it acts indirectly, exchanging energy with the atmosphere constituents visible to the camera FOV. 

Note that the atmosphere, the scattering of the sun on the atmosphere and the sky 
represent the same entity, but here they were separated didactically in order to identify clearly the 
atmospheric irradiance involving a specific path, the source irradiance and the atmospheric 
attenuation. 

In order to see more clearly the interaction of the radiation, it is worthwhile to depict 
what is incident on the target (Figure 3. 4) and on the sea (Figure 3. 5). 

In these levels of presentation, no information was added regarding the dependence of 
radiance on geometric and other variables such as temperature and wavelength. But those 


variables have a great impact on the values of the radiance and should also be considered. 
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Figure 3. 4 Radiation Hitting the Target: Sources and Attenuation. 
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Figure 3. 5 Radiation Hitting the Sea: Sources and Attenuation. 
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C: TEMPERATURE DEFINITION 


Planck’s law relates black body irradiance to temperature and wavelength.: 


G 
M(A)= Se (3. 1) 
eAeT _| 
Where: 
4 
° G37 10 
cm 
° C> = 1.4388 10° ume K 
- wan ee 
cm~ @ Lim 


Then, the apparent temperature of an object (7,) 1s, by definition, the corresponding 
temperature of a black body, that would radiate the same amount of energy as that radiated by the 


object [Ref. 2]. 


[EM (4) bee 7 [EM (A) | Black Body at Temperature Ty, ) (3. 2) 


Applying the Planck law: 
is l 
1 
[EM (A) Laiect = 3° aye (3. 3) 
Al ei we 


For opaque objects the expression for the radiant exitance can be written as the sum of its 
thermally emitted exitance plus the reflected exitance components of radiation (R) from different 


sources, thus: 


A2 





Ci é 
EM(A)Lajece = |] * at R e (3. 4) 
Al A AeT 
e~ “9 —| 

Where: 
e € is the emissivity. 
é To is the object temperature. 
e R is the reflected exitance. 


It should be clear that the radiant exitance of the object will be attenuated by the 
atmosphere and it will combine with other exitances (Figure 3.3) before it reaches the sensor. 
Hence, the final expression will be very complex, and some simplification will be required to 


solve it. 


1. Ship Radiant Exitance 
Considering the ship as the object and considering it as an extended Lambertian source, 
then: 
M=nreL (335) 
This assumption can be made because the ship image appears to have a flat and uniform 


radiance distribution and in consequence, using Figure 3. 4: 
R= melo, 5¢Lo + py5%Ts *Ls + Pr sky *Lsky) (3. 6) 
In the above equation (;4, P;,5 and /; 5%, each represents the reflectivity of the ship 


for the radiation from the sea, sun and sky respectively. It is a function of the direction of 
incidence of radiation on the target and the wavelength. 


Using Figure 3. 5, the sea radiance “Lo” can be expanded one level more. 


Lo = Lsea+ Lrsky + Ls Gry 
Where: 
© Lsea is the sea emission. 
8 Lrsky is the reflected sky radiance. 


Shumaker [Ref. 2] wrote that most military targets approximate gray bodies in the 
8-12 um band with emissivities of the order of 0.9 to 0.95 in that spectral region. The emissivity 
of the R/V POINT SUR has been reported as 0.95 in a private communication (A. W. Cooper), 
and in consequence Is in accordance with Shumaker. 
For an opaque target in thermal equilibrium the emissivity and the reflectivity are related 
by: 
E+p=l1 (3. 8) 
Thus high emissivity means low reflectivity, the value of “R” in Equation 3.4 will be 
small and the major contribution to irradiance from the ship will come from its emission. Hence, 
the target reflectivity will be constant and the maximum error that will result from this 


assumption will have the magnitude of the reflectivity. 


a2 


P= Pt.0 = Fis = Pasny 0.05 (3.9) 


R= ep, elo +7, els + Lop) (3. 10) 


2 Sea Radiant Exitance 
If the object were a source other than the ship, the “R” expression would change. 


Assuming again that the sea is a Lambertian source and using Figure 3. 5: 
R=rel(p,¢Ley +t; eL;) (real) 


For the sea, any energy transmitted through the surface is considered absorbed (except 
perhaps in very shallow waters) and therefore we may write: 

E+ p=\| GF i2) 

Of those backgrounds considered in Figure 3. 3, the sea is the hardest to solve for the 
complicated geometry of reflection due to the randomness of wavelets and its significant 
relationship with atmospheric conditions, mainly the wind direction and intensity. But, it has 
been studied extensively, and the SeaRad program can be used to help in computing the exitance 
contribution from the sea. 

The sea emissivity 1s a function of the wavelength and the radiation incident angle. But 
from [Ref. 2], the use of the average emissivity introduces an error less than 2% for the 4-10 um 
band. 

The sea emissivity does not depend on the sun position but it has a large dependency on 


wind intensity and direction and on the weather conditions. 


=. Sky Vault Radiant Exitance 

From [Ref. 11], the radiance of the sky is the sum of the emission by the atmospheric 
constituents and the scattering of the sun’s radiation. Molecular scattering is important only in 
the short-wave region neighboring the visible, and the emission is important for wavelengths 
longer than 4 um. 

Note that the reflectivity concept is now replaced by the scattering concept, which is 
more convenient to deal with for the aerosol and molecular atmosphere components. SeaRad 
does not make any distinction between the emission and the scattering radiance, and refers to the 


sum of the two radiances as “path radiance”. 
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4. Atmospheric Radiant Exitance 
Here, the same considerations from the previous Sub-subsection 3 will be applied, but 
now, the path considered is the one that contains the air mass between the sensor and the object 


under study instead the one between the sky and the sensor. 


2 Sun Radiant Exitance 

The sun radiation will heat the atmospheric constituents, adding a contribution to the 
atmosphere emission, and also will generate three kind of radiances: the scattered sun radiance, 
the sun radiance specularly reflected from the sea surface, and the sunlight scattered by particles 
beneath the sea surface. The latter, [Ref. 9], 1s nearly independent of the presence of waves and 
can be neglected for large distances from the observer. 

The scattering 1s a complex phenomenon that involves atmospheric constituents, optical 
and solar path, sun position and the model of scattering: single or multiple. In [Ref. 5], the 
authors claim that a single scattering model 1s good only when looking near the sun, since the 
scattering then is dominated by forward peaks. 

The sun radiation reflected from the sea surface is also very complex since one of its 
variables is the wavelet geometry, which is a random variable. 

Fortunately, there is a program named Sea Radiance Prediction Code (SeaRad) that 
manipulates those complexities, and there is no need to worry about the “R” expression for the 


reflected sun radiation on the sea surface. 
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IV. INPUT PREPARATION 


This chapter will present the available data and some data manipulation that will make 
them ready to be used as the input for the range compensation method. 

The core of the work is based on the utilization of SeaRad either to justify a point of view 
or to get new data. 

To obtain the sun coordinates a shareware version of the Cybersky 3.0.2 software was 
used. For mathematical calculations or image processing Mathcad or Matlab in combination with 


the PIRWIN software was used. 
A. BASIC INPUTS 


The basic inputs are those that define the problem. The file “NSD1003.ptw” was chosen 
as a model for calculation of range compensation because it shows the ideal image conditions for 
this analysis. The target is in the center of the picture, the target contour is well defined (good 
contrast), and the horizon line 1s sharply defined. 

The main data list and its origins are listed below. 

From [Ref.8]: 

e Observer latitude: 32° 39’ 55” = 32.6653°. 

e Observer longitude: 117° 14? 31” = 117.2419°9. 

From [Ref. 7]: 

e Observer height: 27.4 meters. 

From the image file “NSD1003.ptw”: 

e Date: April 10, 1996. 

e (Hour:minutes:seconds): 10:27:28 AM PST = 17:27:28 GMT 

2 Aw = (176 — 152) =24 lines (see Figure 2. 15). 

& w =256 lines. 


From the meteorological file “eop30sec.101”: 


’ Wind speed: | m/s. 

5 Wind direction: 238°. 

- Air temperature: 16.4 °C. 

. Relative humidity: 73%. 

e Air pressure: 1,016.2 mbar. 


aS 


e Sea surface temperature: 17.6 °C. 


. Ship latitude: 32° 39’ 6” = 32.6517°. 

2 Ship longitude: 117° 13° 42” = 117.2283°. 
From Cybersky: 

e Sun altitude: 48° 58’ 8” = 48.969°. 

e Sun azimuth: 118° 29’ 54” = 118.498°. 


From calculation, using Equation 4.1 through Equation 4.6: 
° Target distance from observer: 1,980.7 meters. 
e Target Bearing: 139.91° East from North. 


For the determination of the target range and bearing the following equations were used: 


ao aaa (4. 1) 
oe 

b=—- 4.2 
5 Ho (4. 2) 

a (4. 3) 
uaa : 

C= cos'(cosb cose +sinbe since cos A) (4. 4) 

range =reea (4. 5) 

Ling = sini Sessile ve) (4.6) 

sina 


Here re is the radius of earth, the observer position ts (15,0) specified in latitude and 


longitude, the target position is (u,v), also in latitude and longitude, and the bearing is 
measured taking the North Pole direction as the reference position. 
From Equation 2.13 and using the 7° lens, x7 = 1,917 meters. This shows that: 
© FOV used for this image is 7°. 
From experiment notebook: 
e Weather: overcast with cloud and 3,500 feet scatter. Thus, the weather input is 
ICLD = 2: Altostratus cloud with base 2.4 km. 


The geometry of the experiment for this particular case is shown in Figure 4. 1. 
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Figure 4. 1 The Geometry for the File “NSD1003.ptw” Case. 


Note from Figure 4. 1 that the Sun is in a good position for glint formation, but 
unfortunately the wind intensity is very low. Thus, the wind does not favor high values of Sun 
glint. However, the main limitation for setting up the Sun glint is the overcast condition on the 


experiment date. 


B. BASIC CONSIDERATIONS 


I. SeaRad Input File 

SeaRad uses “Tape5” as its input and the format is nearly the same as its predecessors 
LOWTRAN and MODTRAN. The variation is mainly in the additions that SeaRad made in 
computing sea radiance. | 

In this thesis the following basic instructions for SeaRad will be used: 


° Run option: MODTRAN. 


® Atmospheric model: Midlatitude Summer. 
® Atmospheric path: Vertical or slant path to space. 
e Program execution in radiance mode with solar scattered radiation. 


e Program executed with multiple scattering. 


. Use of U.S. Standard altitude profile. 

° Boundary temperature: Sea temperature. 

* Surface albedo (ocean): -5. 

e Type of extinction for boundary layer: Navy Maritime extinction. 

e Seasonal aerosol profile for boundary layer: Spring-Summer. 

° Type of extinction and aerosol profile for stratospheric-layer: Background 


stratospheric. 


e Air mass character: 3. 

e No use of the Army Vertical Structure Algorithm for aerosol boundary layer. 
e Use of default meteorological range. 

s No rain. 

® Observer altitude: 27.4 meters. 

e MIE-generated SeaRad internal database of aerosol phase functions. 

Ze Sea Emission © 


Because the sea is the main background in the EOPACE experiment, most range 
compensation will happen on the sea image. Thus, it is required to find the sea radiance behavior 
as a function of the incident angle (zenith angle of the observer; seen from the footprint). Using 


SeaRad the data of Table 4. 1 can be found for the file “NSD1003.ptw” boundary conditions. 


Incident | Sea Emission Average |Sky Reflection| Sun Glint Total 
Sea (watts/m7/sr) | watts/m*/sr irradiance 
Emissivity ae IS) 

89.8 | 14.1688 | 0.4234 | 19.6022 = ee 
. 19.5307 ee ens 
aC 
a 
rs | ie.si40 | 0se5_ [a7 [0 —«* 082i 
2 2 
[a0 [33.6995 [o.70s0[93sser 033.0502 
7s _|_273tst_| o8isa [sara [0 —~«Y—8.2a892_TS 
roa 8695 0.8850] seas | 0 3908 | Ts 
65s [_31.05535_| 0.9304 | 22sa21_[ 0 | 3328956 [<7 
60 | si8s2it_[ 0.956 [143539 | _0 | ~33.28750_ | 17S 
(30 | _32.58503_| 097e | oa | _0 | 3329877 [__175__ 


Table 4. 1 Sea Radiance Data for File “NSD1003.ptw”. 
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To get a feeling for the influence of weather on the reflectivity, that 1s (1 - emissivity), 


and to compare with other authors’ curves, the plot of Figure 4. 2 was built. Unfortunately the 


authors of [Ref.2] and [Ref. 6] did not specify the weather conditions. However, from 


Figure 4. 2 and considering the magnitude of the results, it is clear that both authors considered 


the “no clouds” condition. 


Notice the effect of wind on the reflectivity: as the reflectivity increases the wind 


intensity also increases. 
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Figure 4. 2 Comparison of Reflectivity. 


Using polynomial fitting, the sea reflectivity “po, ” for 1 meter/sec wind speed of this file 


can be written as a function of the angle of incidence “a” in degrees as: 


p,(a)= D6ea® + D5ea?+D4ea* + D30a>+D20ea*+Diea + DO (4. 7) 


Where: 
eo 


Do: -74.9072 
DI: Gia 


a2 


e Dz: -0.269547 


© D3:  0.00553385 

e D4: -6.36474E-005 
e DS: 3.88588E-007 
. D6: -9.81879E-010 


This formulation will be used in Chapter 5 for correction of the SeaRad sky reflection. 


a Picture Geometry 

One important factor in atmospheric attenuation is the distance. Thus, it is necessary to 
determine all the points in the picture that have equal distances relative to the observer. 

Considering a plane earth with plane equation: z = 0, and the observer positioned at the 
rectangular coordinate (0, 0, h), then the locus of all points on the earth equidistant from the 
observer will be the intersection of the sphere of radius re (all points in the space that have the 


same distance to the observer) and the earth, which can be written: 


x? ty? = re* —h* (4. 8) 


This shows that the observer will perceive the points having the same distance from him as lying 
onacircle. The geometry is consistent with Figure 2. 14. 

Assuming that the nearest target picture will be at 1 km and that the target is in center of 
the picture, then the nearest point in the picture will be at approximately 308.83 meters from the 
observer (see Figure 4. 3. b). Even for this small distance the points at the same distance to the 


observer can be approximated as lying on a horizontal line with at worst an error of 0.58 meters. 
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308.83 
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a) Elevation View a) Plan View. 3.5° is the azimuthal 
angular width of the scene 


Figure 4. 3 Picture Geometry: General Case. 


This fact is very favorable for the solution of this thesis since in order to make a picture 
compensation all that is necessary is to know the atmospheric parameters as a function of distance 
for just one particular azimuth direction, specifically, the target direction. 

The ship itself should be considered differently from the sea because all points located on 
the ship can be considered to have the same distance to the observer. That fact will be beneficial 
because the ship contrast will be improved when using the range compensation method where the 
sea surface background approaches the horizon (the background temperature will decrease). 

For figuring out the incident angles of interest for the file “NSD003.ptw”, Equation 2.10 
can be used to build up Figure 4. 4. 

The target-horizon distance in image coordinates is: Aw = (176 — 152) =24 pixels. The 
distance from the lowest point on the image plane to the horizon is: Aw = (256 — 152) = 104 
pixels. The target position considered was location 3 of Figure 2. 16, that is “water aft” probe 


position. 
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Figure 4. 4 Picture Geometry in Vertical Plane for File “NSD1003.ptw”. 


Considering that the zenith angle of the observer line-of-sight is 180° minus the incident 
angle, for the sea compensation this angle will vary from 93.02° to 90.168° and for sky from 
90.168° to 86.02°. The horizon zenith angle of 90.168° is only an approximate value. The 
correct horizon zenith angle can be determined by plotting the sky path radiance versus observer 


zenith angle, as will be done later in Subsection 5. 


4. Target Temperature 

The target itself is the main component in the thermal image picture. Therefore, it is 
desirable to check the quality of the ship temperature distribution recorded with the AGA, Figure 
4.5, by comparing it with the skin temperature recorded with the thermistors, shown in 
Figure 4. 6. | 

In fact, it is impossible to associate quantitatively each channel position with its 


corresponding point in the picture due to the following facts: 


e The ship is not perfectly transverse to the target line-of-sight. 
° The position coordinates of the thermistor are very loosely specified. 
® The exchange of thermal energy between the ship and the surroundings produces 


temperature perturbations around the ship. Thus the exact ship profile cannot be 
determined from the image. 
Nevertheless, qualitatively the data seem to be correct, with a region of high temperature 


near the port stack, another region with low temperature that follows the ship profile and finally 
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the surrounding sea temperature. The comparison should also be considered qualitative because 
the AGA temperature is an apparent temperature that incorporates the path radiance, atmospheric 
radiance, radiance reflected by the target, and the target emission radiance. On the other hand, 
the thermistor temperature considers only the target emission temperature. 

Table 4. 2 can be built for these three regions, where the 2.3+ 1.7 °C added to the AGA 


temperature refers to the calibration correction. 





Region 








Table 4. 2 Qualitative Comparison of Temperature: Corrected Radiometric Temperature Versus 
Thermistor Temperature. 





Figure 4. 5 Thermal Image Temperature Distribution Uncorrected for Path Radiance. 
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Skin Temperature: April 6, 1996 
10:27:28 AM PST 


11-16.386 8-16.754 14-26.857 4-16.401 


10-17.267 216.918 


3-14.572 


Figure 4. 6 Thermistor Temperature: Temperatures Recorded by Selected Thermistors in File 
“eop30sec.101”. 


Note that in file “eop30sec.101” the METOC measurement of sea surface temperature 1s 
17.6 °C, whereas the skin thermistor number 3 recorded a sea surface temperature of 14.5 °C. 
Comparing these outputs as time series as shown 1n Figure 4. 7, it is clear that the METOC value 


follows the thermistor value with a time delay. Several reasons can be advanced for such 


differences. 
The influence of this difference in the sea temperature will be analyzed better later in the 


following Range Compensation Chapter. 
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Figure 4. 7 Comparing Temperature Data: Sea Surface Temperature from Thermistor 3 
Compared with R/V POINT SUR METOC Measurement. 


= The Horizon Determination 

The correct identification of the horizon position is needed since the sky path radiance 
has a pattern completely different from the sea path radiance. Figure 4. 8 shows computed path 
radiance versus frequency with zenith angle as a parameter. We see that at any frequency the 
path radiance rises as the horizon is approached from either direction, with a sharp maximum at 
the horizon. This condition is represented by the dotted line for 90.15°, identifying the optical 


horizon. 
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Figure 4. 8 Path Radiance for Different Zenith Angles. The Maximum Path Radiance at 90.15° is 
Indicative of the Optical Horizon. 


Running SeaRad for other observer zenith angles and integrating the path radiance versus 


frequency curve, Table 4. 3 can be built for varying zenith angles. 


Zenith Angle | Path Radiance 
78.25 33.4 


90.16 Sony 
- 90.20 B0R7 


90.60 
50.80 


Table 4. 3 Spectrally Integrated Path Radiance as a function of Zenith Angle 
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Plotting the data of Table 4. 3 will provide a better view of the horizon position as long as 


the path radiance of the sea is much smaller than the sky path radiance (see Figure 4. 9). 
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Figure 4. 9 The Horizon Determination. The Sharp Discontinuity in Radiance Identifies the 
Horizon. 


Finally, from Figure 4. 9, the observer will see the horizon at a zenith angle of 90.15°. 
This horizon zenith angle takes into account the ray bending and the spherical earth. Using plane 
earth calculations 90.168° was found instead of 90.15° (see Figure 4. 4). This shows the validity 
of using plane earth geometry. 

In terms of temperature, the sky path radiance can be seen better in the isotherm display 
mode as in Figure 4. 10, where the calibration correction has not been applied yet. Note that 
unlike the sea the sky does not have random sources of radiance such as the ones seen in the 


lower left of the image, unless broken clouds and sun glitter are present. 
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Figure 4. 10 Uncorrected Sky Isotherms. 


According to MacCartney, [Ref. 3], the calculation of the sky radiance imposes two 
difficulties: the sky does not have area and, principally, the sky does not have definite distance 
from the observer. This means that the concept of range compensation loses its meaning when it 
is applied to the sky because the distance is indefinite. 

For that reason, the range compensation will apply only for the sea and the ship with the 
observer zenith angle varying from 93.02° to 90.15° for this particular file “NSD1003.ptw”. 
However, the sky radiance will not be overlooked and will play a major role in tuning up the 
SeaRad sky radiance by adjusting the sky path radiance and comparing the SeaRad output with 
the temperature of Figure 4. 10. 

To determine the sensitivity of path radiance to the major parameters used in SeaRad, the 
comparative Table 4. 4 was.constructed from multiple SeaRad runs with varied values of Air — 
Mass Character (ICSTL), Cloud and Rain Model (ICLD) and Surface Meteorological Range. 
The calculated results show that changes in ICSTL and in Surface Meteorological Range produce 


insignificant variation in the path radiance over the range of parameters used. 
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Table 4. 4 Influence of Meteorological Parameters on Path Radiance for the Horizon Zenith 
Angle. 
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From Table 4. 4, it can be seen that there is a lower limit for the path radiance below 


which the sky path radiance cannot be decreased. 
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V. RANGE COMPENSATION 


This chapter will present a step by step procedure for the range compensation method for 
the file “NSD1003.ptw’. As stated in the previous chapter, the range compensation will be 
applied only for the sea and the ship, with the observer zenith angle varying from 93.02° to 
Dalam 


AS DATA MANIPULATION 


The range compensation process is not an easy task and the complexity can be 
appreciated with the help of the picture of Figure 5. 1. For each specific distance the observer 
measures an apparent black body temperature of Tapp-FR. This temperature corresponds to an 
irradiance Ms (area under the black body irradiance curve) that allows the paleelation of the zero 
range irradiance Mo, knowing the atmospheric transmittance and the path radiance. The range 
compensation then ends for that distance by computing the apparent black body temperature at 


zero range, Zapp-ZR, which corresponds to Mo. 
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Figure 5. | The Range Compensation Problem. 


>i 


Note that the above process is very extensive because it should be applied over all pixels 
of the image plane and in consequence a better algorithm should be found. 
The first step in that direction is to define the tools that will be used: matrix variables. In 


order to manipulate the wavelength dependency for ‘the infrared radiation parameters, first, the 
smallest frequency step size allowed by MODTRAN (5 cm! ) needs to be selected. Then, all 
SeaRad output should be standardized to vary over the same frequency bandwidth of 830 cm! 


to 1250 cm! , with a total of 85 frequency points in this range. That is: 





830 yi 
835 i) 
eal eee: (Say) 
1245 Vg4 
[1250] | vgs | 
12.048 Ay 
lle /6 Ay 
~- | 
= eo) = =e (5. 2) 
V 
8.032 Aga 
| 8.000 | | Ags 
Where: 
« V is the constant frequency vector in cm. 
° A is the wavelength vector in zon. 
10,000 <a 
e — represents a Matlab operation of scalar division by elements of the 
V 
frequency vector. 


The result of this standardization is that the SeaRad output can now be manipulated as a 


matrix, as shown in Figure 5. 2, where each of its columns shows a frequency dependence. 


a2 


SeaRad Output File (SOF) Matrix 


8.32E-01 830 0.346 8.32E-01 
1.39E+00 835 0.345 1.39E+00 
EY3E700 840 0.347 1.93E+00 
2.48E+00 845 0.346 2.48E+00 

* * * * 

* SOF = * * * 

* oe » * * * 
3.96E+01 1235 U. 233 3.96E+01 
3.99E+01 1240 0.22 3.99E +0!) 
4.02E+01 1245 0.202 4.02E+01 
4.04E+01 1250 0.18 4.04E+01 





Figure 5. 2 Transformation of Text File into Matrix Variable. 


For example, if the third column represents the path radiance, then it can be written: 


Lpath(v,)| [8.32E-01 
Lpath(v>)| |1.39E +00 

Lpath(v) = SOF (*,3) = = (5. 3) 
Lpath(vg,)| | 4.02E+01 | 
Lpath(vgs)| | 4.04E+01| 


Accordingly, all data manipulations will be treated as matrix operations. 


Frequently the average value is needed. Its integral form is as follows: 





b 
average of flv) = 7 ® |#v)av | (5. 4) 


i= 
Taking again the path radiance as an example, we will have for the average path radiance 


in the interval (830 cm: 1250 cm!) in matrix form as follows: 


Av; 


Av 
[Lpath(v;) Lpath(v>) -- Lpath(vgq) Lpath(vas)J] <° | (5, 5) 


(Gent. = 
ae (1250 — 830) 


Avs4 
LAvgs 


Another term that it is commonly used 1s the Integrated Radiance, defined as: 


aS 


Av; 
Av» 
Lpathro; = [Lpath(v,) Lpath(v>) -- Lpath(vs x) Lpath(vg5 eb ; (S. 6) 
Av34 
[Aves | 
Having found the total radiance, Equation 3.3 can be used to find the apparent 
temperature. The solution of this equation should be done interactively as was done in subroutine 
RtoT() of SeaRad. 
Subroutine Rto7( ) in SeaRad computes the apparent black body temperature 
corresponding to a selected radiance by an iterative process. Starting with a “guessed” 
temperature, the corresponding output radiance is calculated, and the input corrected until the 


output matches the required radiance within a selected error range. 
B. SKY RADIANCE ADJUSTMENT 


The sky path radiance dependency on the zenith angle alone will impose the condition 
that the sky temperature will also be dependent only on the zenith angle. This means that at the 
horizon, as an example, a unique value of an apparent black body temperature can be found and 
used to compare with the temperature read directly from the image plane and also used to adjust 
the sky radiance. 

Using the default meteorological range, with sea temperature of 17.5 °C, air mass 
character 3, no clouds or rain condition and using the observer zenith angle to the horizon, that ts: 


90.15°, then the execution of SeaRad will result in the plot of Figure 5. 3. 
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Figure 5. 3 Smaller Sky Path Radiance Condition. 


Note from Table 4. 4, that the ICSTL=3 case gives close to the lowest radiance values. 
Using that radiance as input for Equation 3.3, and solving it for the apparent black body 
temperature, a value of 21.4 °C is found. 

Thus, the SeaRad result is higher than the expected values of 18.34 1.7 °C (Figure 4. 10) 
by a nominal value of 3.1 °C. 

Since just changing the SeaRad parameters for this specific case cannot decrease the path 
radiance further, the solution is to offset the path radiance by a constant radiance value. 

Taking the path radiance vector “Lpath” from Figure 5. 3, the required offset can be 


estimated by the following process: 


e Guess an offset. 

° Change vector Lpath to Lpath = Lpath — offset*n; where n is an integer. 

e Apply Equation 3.3 to find the average black body temperature over the 
frequency band. 


SS 


Using SeaRad, the results of Table 5. 1 can be obtained. By interpolation from this table 


the offset needed to correspond to the required temperature of 18.3 °C is found to be 5.76 107 
































watts/cm’. 
Offset Average Temperature 
(watts/cm’) (°C) 
2.2601e-004 2021 
3.390 1e-004 19.63 
4.5201e-004 18.99 
5.6502e-004 18.36 
6.7802e-004 17.78 
Table 5. 1 Finding the Correction for SeaRad Path Radiance. 
This offset can also be written in me using the following equation: 
m” .sr.cm 
path offset = 5.76 1074 — Ss ger agen ieee 7 ———— (6.7) 
cm S5e.zel0 © m~esrecm 


Where: It was assumed Av = Sem~’. Thus, the offset of 5.76e-4 watts/cm” should be 


added to all sky path radiance in all the SeaRad calculations that will follow this section. 
G. SEA RADIANCE ADJUSTMENT 


Like the sky path radiance, the sea radiance can also be adjusted since we usually know 
the sea temperature. The idea here ts to calculate the reflected and emitted radiances that leave 
from a given sea position at a given apparent black body temperature read from the 
meteorological data. Adding the radiances introduced by the atmospheric propagation, the 
apparent temperature at the observer position can be found and compared with the value read 
directly from the image plane. If there is no difference in the compared temperatures, then no 
correction is needed; otherwise a similar procedure to that used in the previous section should be 
applied. 

From Figure 4. 4, the target position relative to the observer has a zenith angle of 90.83°. 
Using this angle and all other inputs from the previous section, we can build the following 85 

watts 


component radiance vectors by running SeaRad, all in 


m~ esr ecm! 
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2 Footprint to observer path radiance: Lp = SUMT. 


° Sun reflected on the sea at zero range (leaving footprint): Lsuno = Nvo. 

e Sun reflected on the sea at full range (reaching observer): Lsun = Nvo.Tau(]). 
e Sea emission at zero range: Lseao = e.Nbb. 

e Sea emission at full range: Lsea = e.Nbb. Tau(I). 

“ Sky reflection at zero range: Lrskyo = Nvsky. 

® Sky reflection at full range: Lrsky = Nvsky.Tau(1). 


In the above radiances, the left side of the expression is the variable definition, and the 
right side is its SeaRad interpretation, where: 

® e = sea emissivity. 

2 Tau(I) = atmospheric transmissivity. 


e Nbb = Planck’s black body radiance. 


Comparing those variables with Equation 3.5 and Equation 3.10, we identify: 





Lrskyo = Po ® Ly (5. 8) 
Lsuno =T,¢L, (5. 9) 
Lseao = aA o —_*sea__ 5. 10) 
Te A> _o2 
ot *Tsea = 


Writing for zero range sea radiance: 
Lo = Lseao + Lrskyo + Lsuno (5.11) 
Considering the sky path radiance correction of Equation 5.7, we will have after its 


reflection on the sea, using the polynomial form of Equation 4.7 for p, : 


(5. 12) 





sky _correc = Po F - 90.83° Je 5.76 @ 
5°77 


Lo = Lseao + (Lrskyo _ sky _correc)+ Lsuno (oe 15) 
And for full range: 
Lo@camera = Loet+Lp = Lsea + Lrsky — sky _correc + Lsun + Lp (5. 14) 
For this experiment date the sky was overcast and in consequence the Lsun component can be 
neglected, simplifying the previous equations. 
As was done in the previous section, by application of Equation 3.3, using Equation 5.14 


as input, the apparent black body temperature at the observer position can be found. The sea 


oT 


footprint position was kept constant at a zenith angle of 90.83° relative to the observer. The 
results are shown in Table 5. 2. 
Sea Temperature (°C) ' Calculated Apparent 


(Tbound) 












Temperature Seen by AGA 

(°C) 
12.6 
17.6 (Metoc) 13.4 


OS 


Table 5. 2 Temperature Attenuation for the Ship Position. 













14.5 (thermistor) 





Comparing the Table 5. 2 results with Figure 4. 5, and considering that the average 
Tbound temperature provides the correct apparent temperature for the sea at the ship position, no 
extra radiance compensation Is necessary. 

To see the influence of the path radiance as the observer zenith angle approaches the 
horizon, Figure 5. 4, for an observer zenith angle of 90.2°, and Figure 5. 5, for an observer zenith 
angle of 90.85°, were plotted. In both figures the same sea temperature Tbound=289.6 TC was 
used. It can be inferred that the increase of path radiance, as the observer zenith angle approaches 
the horizon, is the main reason for the increase in the apparent temperature at the observer 


position. 


58 


Radiance in watts/cm 21m 


2.5 


2 
oe ae Sea 
5 Reflected sky 
Path 
1 
Os 





PO eee 


800 850 900 950 1000 1050 1100 1150 1200 
Frequency in cm-1 


1250 


Figure 5. 4 Radiance Near the Horizon. Observer Zenith Angle of 90.2°. 
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Figure 5. 5 Radiance for Near the Target Position. Observer Zenith Angle of 90.85°. 


THE TEMPERATURE DISTRIBUTION IN THE IMAGE PLANE 


Assuming a plane earth, Figure 5. 6 can be used to find the relation between each 
observer zenith angle a@ and the horizontal projection distance x, or between q@ and the row 
numbers p in the image plane. Here the subscript H refers to the horizon and the subscript N to 


the nearest point that the AGA sensor can see (bottom of the image). 


a7 





Figure 5. 6 Correspondence Between Zenith Angle and Horizontal Distance. Geometry in 
Vertical Plane. 


The data (a N:@,@ H) can be associated linearly with the counterparts (p Nae He in 


the image plane. The p value is an image frame characteristic. The PTRWIN files have 256 lines 
and 125 or 250 pixels per line. Then, p varies from 1 to 256 (1 at the top and 256 at the bottom). 


AH -A2n PH-PN 
And finally write: 
aA-Aan je - 


CoH a GuN. 
Comparing Figure 5. 6and Figure 4. 4 and analyzing the image plane using the PTIRWIN 


software for example, the following data can be found: 


e Pn = 256. 
® Pu = 152. 
. ay = 93.024° 


e QH = 90.168? 


° 90.168° < a < 93.024° 


For each p value there is a corresponding pixel value in digital level. Using Equation 2.3 


each digital level can be transformed into Isothermal Units, and finally, by using Equation 2.1 
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the temperature associated with each p value, that is 7(p), can be found. The plot of 7/p) as a 
function of the observer zenith angle @ is presented in Figure 5. 7, for column 64. Note that the 
calibration correction is not applied yet, and that column 64 corresponds to the direction that 
contains the hottest ship point. There is, the ship hottest point is contained on the vertical plane, 
depicted on Figure 5. 7. 

24.00 — 


| Ship Hottest Point ("Hot Spot") 


20:00" 


Apparent Temperature (Degrees Celsius) 
| 


90.00 91.00 92.00 93.00 94.00 
Observer Azimuth Angle (Degrees) 


Figure 5. 7 Apparent Temperature Distribution Versus Observer Zenith Angle for the Vertical 
Plane (Column 64) That Contains the “Hot Spot”. 


E. DETERMINATION OF THE “REAL” SEA TEMPERATURE 


When we look at the sea near the horizon, the sky reflection dominates the sea emission 
component and the sea temperature will look lower than it really is. That temperature is referred 
to as the apparent black body temperature and it corresponds to the sea “Tbound” SeaRad 
variable. Actually it is an apparent temperature because it includes all radiation leaving the sea 
and not just the sea emission, this latter being associated with the real temperature of the sea. 

To find the apparent black body sea temperature, SeaRad can be used with a little 


modification on the subroutine RtoT( ) to consider the sky correction, Equation 5.12. A 


6] 


dimensional change is necessary, and it can be done by just subtracting the sky correction from 
watts 


2 


the total radiance in . 
cm~ @ sr 


Then, for Rto7() the sky correction 1s: 
l 
sky _correc = (l—esea)e5.76 «— (5. 17) 
1 


Where: esea is the average sea emissivity also calculated from SeaRad, applying 
Equation 5.4 to the emissivity variable “e”. Now, keeping the sea temperature constant and 
varying the observer zenith angle, Figure 5. 8 can be constructed by running SeaRad with the sky 
correction. 

By plotting together Figure 5. 7, corrected for AGA calibration, and Figure 5. 8, and 
using the same coordinates, the apparent sea black body temperature can be read directly from 


Figure 5. 9, from the intersections of the curves. 
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Figure 5. 8 Sea Apparent Temperature at Observer Plane Versus Observer Zenith Angle with Sea 
Apparent Temperature at Sea Plane as Parameter. 
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Figure 5. 9 Finding the “Real” Sea Temperature. 


The above method will become tedious when processing an image with 256 by 125 


pixels, and an automatization is required. 


The idea is very simple: 


Use Equation 2.3 to transform the original image in digital level to isotherm 
units. | 

Use Equation 2.1 to transform the isothermal unit image into a temperature 
image, called a JT image. This image file will have dimension 256 rows by 125 
columns. 

Correct 7 image in order to take into account the calibration error by adding to 
all its elements the nominal error 2.3 °C. 


Define the vector @. Its elements should be in the range 


90.16° <@ <93.02°. Also, find the corresponding p vector, using 


Equation 5.16. 
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90.16° 152 
90.46° 163 
90.76° 174 
91.06° 184 
91.36° 195 
@=191.66°| % p=|206 
91.969 217 
91.26° 228 
92.56° 239 
92.86° 250 
93.02° 256 
e Define the vector of sea Tbound values. Here 11 values were selected, but 


Tbound, p or acan have any size and any values of interest. 


Thound =[- 263.15 -70 -20 0 12 14 16.5 19 21.5 24 40] 
. Now, as before, by running SeaRad with the sky correction for the sea 


temperature Tbound and the observer zenith angle a, Figure 5. 8 can be plotted 


again or the following matrix built: 


17. @ 7.9 mee. 4 1047 lo 16,9) 19° {oi 1/9 i>) 192 tee 

ae 4.7 9.1 ee 14 4g (OID CIS 623 ee 
— I ee 2 es 122 12.7 19.3 14 14:7 15.4 ae 
=—/.5 =<$=49 29 78S 11.2 118 12.6. 13.3 14d- 449g 
—10.5. =7:7 welezo 7 0 Sia 1226 S./ 146 
—15) 9m Uf Jomo we 6 6II IZ 13:4 14.4 2 
—o4  — 12 1.4 D4” OZ Nall 12.1 igh 4] eee 
7 3 eee WO. 11 PRE TS 2. 14.2 274 
= [9 tems. 3 a= 30 8] eS. Same iam 10.5 ae mel 2.6 «13.8 212 
—21.5 -17.4 -48 2.7 7.7 86 97 108 21.9 13.1 20.8 
Sl EO ae ees ie OZ Vial oe! 6 12.8 20m 


Each column “J” of Tapp(*,J) corresponds to Tbound(J), and each row 


Tapp(I, *) corresponds to p(J) or equivalently to a(Z). 


Where: * means all elements. 


=e of 


The key point here is that at each p(I) coordinate, Tapp/(I,*) can be applied over 


all pixels in all columns of this row (see Chapter IV.B.3). Then the problem 


resumes in finding the vector Tapp_row by interpolation, corresponding to each 
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integer row number p’ of the image plane from 152 to 265. As an example, for 
row p =152: Tapp _row = Tapp(1,*), and for row p =153: Tapp _row should be 
interpolated between Japp(/,*) and Tapp(2,*), knowing that in this interval p(/) 
varies between 222 and 152. 

+ The Tapp_row value obtained for the specific row p’ can be applied for each 

T image(p’,*) element to find the corresponding Tbound element. This is an 
easy task because each column of Japp is associated with a Tbound element. 
As an example, assume p =152. Then Tapp_row = Tapp(1,*), and interpolation 
over all 125 elements of row 152 of TY image matrix is needed, that is, 
T_image(p’,*). To interpolate a hypothetical T_image(p’,20)=18.8 TC, for 
instance, it can be seen that 18.8 is between Tapp(/,4)=18.7 and Tapp(1,5)=18.9 
and the corresponding Tbound varies from T7bound(4)=0 to Tbound(5)=12. 
Thus, the corresponding sea Tbound value for the apparent temperature of 
18.8 °C is 6°C. 

In the above procedure, from [Ref. 2], the sky was considered as having an apparent 
black body temperature of -40 °C. This is necessary to put a lower limit in the calculations. 

Also, the range compensation for the ship pixels should be discarded because their 
distance to the observer is almost constant, and the ship emits a different Spectrum of radiation 
from the sea. The ship will be range-compensated in the next section. 

The ship temperature is very prominent in Figure 4. 5 with respect to its surroundings. 


Thus, the ship pixel coordinates can be found easily and a mask can be built as in Figure 5. 10 





Figure 5. 10 Ship Mask. 
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The final range compensation, missing only the ship correction, is shown in Figure 5. 11. 
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Figure 5. 1] Range Compensation for the Sea. 


Another ship can be seen on the right side of Figure 5. 11. Its range compensation is also 
wrong for the same reason as for the main target discussed previously. But, it does not have any 


kind of relation with this thesis, and in consequence it will not be corrected. 
F. RANGE COMPENSATION FOR THE SHIP 


From Figure 3. 4 and Equation 3.10, Equation 5.18 can be written: 
L, =target emission + p, @ (Lo Sayan Tie Ls) (5. 18) 
Ls is almost zero because of the overcast condition, and before any radiance correction 


Lo is given by Equation 5.11 and Lsky by Equation 5.8. Hence, calling Lte the target_emission: 
| L 
Lig = Lteo + p; (Lseao BLT SKV Oca Lsuno )+ — (5. 19) 
—e 


Finally, considering the sky correction from Equation 5.12: 


Lrskyo — sky _ correc 


It _reflo=p,e (Lseao + Lrskyo — sky _correc + Lsuno )+ ; (5. 20) 
—e 
Lio = Lteo+ Lt _reflo (5.2) 
L,@camera =L,,¢t+Lp = Lie+Lt_refl+Lp (5. 22) 
7 = 
Lt_refl=p,e (Lsea + Lrsky — sky _ correc + Lsun)+ it (5. 23) 


l-e 


66 


All the components of Lt refl, Lp and t can be obtained using SeaRad, as was done 
before. For the target emission calculation the ship average temperature can be found using the 
ship mask. 

But, a better solution is to find the radiance Li@camera corresponding to different 
temperatures on the target and apply Equation 3.3 to find the apparent target temperature at the 
observer position. Comparing with the previous sea procedure, this 1s equivalent to finding just a 
one row Japp vector for the target. 

To find Lt_refl, Lp and t, SeaRad must be run just one time because all the ship points 
are assumed to have the same distance to the observer, and Lt _refl, Lp and t, are functions of the 
sea temperature instead of the target temperature. The main input parameters for SeaRad are: 

® Tbound = 16.04 °C. 

a Observer zenith angle = 90.83°. 


Now, applying a target version of Equation 5.10 to find the target emission, for 


Eship = 9.95 , from Chapter III.C.1. 





Lico= i E a (5. 24) 
ren 7 2 
e °T ship a ] 


For a target temperature of 12 °C, the radiance distribution in frequency for the ship is plotted in 


higuresswl2. 
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Figure 5. 12 Observed Radiance From the Ship Position. The Reflected Components are Small 
Compared to the Emission, but the Contribution of Path Radiance is Significant. 
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Applying Equation 3.3 to Equation 5.23, for different ship temperatures, the 


corresponding in-band average apparent temperature at the observer position can be found. The 
results are in Table 5. 3. 






Target Temperature 






Target Apparent 








Tbound (°C) Temperature (°C) 


Table 5. 3 Mapping of Apparent Ship Temperature to Ship 7bound Temperature. 
Applying the ship mask, 
Figure 5. 10, to the AGA image, only the ship pixels can be extracted. Applying Table 5. 


3 to those pixels, the mapping of the apparent temperature to the target 7bound temperature can 
be done. 


Figure 5. 13 shows part of the data of column 64 of the compensated image, after the 
mapping of the ship apparent temperature to 7bound temperature. 
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Figure 5. 13 Comparison of the Range-compensated Temperature and Original Apparent 
Temperature Data for the Ship before any Sea-Ship Interface Correction. Image Plane Sectioned 
at Column 64. 


The superposition of the range compensation for the ship and for the sea , for the same 


column 64 of the compensated image is represented on Figure 5. 14. 
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Figure 5. 14 The Discontinuity of Temperature at Sea-Ship Interface. Image Plane Sectioned at 
Column 64. 
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Figure 5. 14 shows that an adjustment is needed to correct the discontinuity in the 
sea-ship interface. This continuity should occur only on the right interface side, that is, the region 
nearest to the observer, where the observer sees the hull of the ship encounter the sea. 

The procedure to adjust the ship temperature ‘is: 

e Define a column that contains ship data. As an example, take column 64. 

2 In that column, find the rightmost row number, correction_row, where the ship 

temperature is bigger than zero. 

* The temperature correction is: [Sea temperature at row (correction_row+J])] — 

[ship temperature at row correction_row}. 


The result is shown in Figure 5. 15 and Figure 5. 16. 
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Figure 5. 15 Ship Temperature distribution after Discontinuity Correction at the Sea-Ship 
Interface. Image Plane Sectioned at Column 64. 
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Figure 5. 16 The Final Image After the Application of the Range Compensation Method. This 
Image represents the Apparent Black Body Temperature at Source Origin Position. 


G. ANALYZING THE RESULTS 


A visual comparison of the compensated and the original image, Figure 5. 17, shows that 
a better way to analyze the images must be found. Here, it was chosen to section the image in the 


vertical and horizontal directions. 





Figure 5. 17 Left is the Compensated Image, Right is the Original Image. 


In the vertical direction, Figure 5. 18, the ship temperature has become more flattened. 


This is due to the fact that for small target Tbound the path radiance has a great influence on 
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temperature attenuation and it is responsible for the small values of the temperature differenceAT 
between Tbound and its corresponding apparent temperature. Since the atmospheric parameters 
are kept constant for the ship compensation, and the only changing factor is the ship emissivity, 
then as Tbound increases its influence becomes dominant over the path radiance. In consequence 
there will be a temperature at which there will be a transition from one situation to another and 
where AJ = 0. From Table 5. 3, that temperature is near 23 °C. Since most of ship data are 


around 23 °C this explains the flatness of the ship compensation results. 
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Figure 5. 18 Comparing The Results of The Image After the Range Compensation. The Ship 
Temperature Distribution Has Become Flattened and the Gradient of the Sea Temperature Near 
Horizon Very Accentuated. 


Also note from Figure 5. 18, how accentuated is the path radiance on the sea AT near the 
horizon, where the sea temperature approaches the sky temperature. 

Far from the horizon and near to the observer the path radiance is smaller, so that less AT 
correction is necessary. 

In image row coordinates, the ship is near the horizon and because AZ there is higher, it is 
expected that the contrast in the horizontal direction will be improved. This improvement can be 
seen by sectioning Figure 5. 16 horizontally at different rows containing the ship pixels, as is 


done in Figure 5. 19 to Figure 5. 22. 
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Figure 5. 19 Ship-Sea Contrast for Row 172 (Position Far From the Horizon Row: 152). No 
Visible Improvement in the Contrast. 
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Figure 5. 20 Ship-Sea Contrast for Row 169. Still No Visible Improvement in the Contrast. 
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Figure 5. 21 Ship-Sea Contrast for Row 163. The Contrast is Improved as the Ship Pixels are 
Contrasted with the Sea Pixels Near the Horizon. 
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Figure 5.22 Ship-Sea Contrast for Row 159, a Position Very Near to the Horizon (Row 152). 
The Improvement in the Contrast is Intensified. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The equipment calibration is the first difficulty in working with thermal vision data. The 
correct calibration will allow the correct transformation from the equipment-dependent 
parameters to the real world temperature scale. But in order to find that calibration curve a 
calibration data point associated with the experiment data should exist. That is, black body 
measurements interpolated between the image measurements should exist. Having those data to 
check, it should be easy to adjust the controller output to reset it to the conditions of the time of 
the experiment and find in the laboratory the correct calibration curve. Then, using this correct 
calibration and using the same black body measurements the error and the dispersion associated 
with the experiment can be found. 

In this thesis, an error of -2.3 + 1.7 °C was found for unpolarized images. This 
magnitude represents about a 10 % error in the mean. This error was used to offset all the image 
pixels after the application of the correct calibration curve. 

To develop a range compensation method is not an easy task due to the complex 
mechanisms of the radiation interaction, the randomness of sea wavelets, and because it requires 
knowledge of the theory that involves atmospheric constituents, optical and solar path, sun 
position and the model of scattering. Then, the second step in the solution of the problem 
involves same simplifying assumptions. 

The following simplifications were assumed: 

° All sources were extended Lambertian sources. This allows introducing a 


dependency on solid angle by just dividing or multiplying by 7. 


‘ The emissivity of the research ship R/V POINT SUR was considered constant 
and equal to 0.95. 
® The reflections of the sun radiation from the sea and from the target were 


neglected due to the overcast condition of the experiment. 

° To skip the hard mathematics behind the effects of atmosphere on the generation 
and attenuation of radiation, the SeaRad program was used. 

e Due to the small field of view of the equipment and the big distances involved, 
the geometry of the problem allows considering all the horizontal pixels in the 


image plane to have the same distance to the observer. Then, it is necessary to 
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find the atmospheric parameters as a function of distance in just one particular 
azimuth direction. 

After these simplifications the structuring of the range compensation method was begun 
by considering that all points located on the ship have the same distance to the observer and that 
the sky has an indefinite distance to the observer so that the concept of range compensation has 
lost its meaning. But this same property that the sky radiance is independent of distance allows 
correction of the SeaRad sky radiance by comparing it to the apparent temperature in the image. 
For the file “NSD1003.ptw”, the offset needed to correct the SeaRad sky radiance output is 5.76e- 
4 watts/cm’. 

The file NSD1003.ptw was chosen as the example for calculation of range compensation 
because it shows the ideal image conditions for this analysis. That is, that the target is the center 
of the picture, that the target contour is well defined (good contrast), and the horizon line is 
sharply defined. 

Using SeaRad to calculate the apparent sea temperature Japp _sea at the ship position and 
comparing it with its corresponding value read from the image file, it was found that no extra 
radiance correction was necessary for the sea radiance. 

The sea can then be range-compensated by first building curves of apparent temperature 
Tapp_sea versus observer zenith angle a@ for different sea apparent black body temperature 
Tbound_ sea. Next, by finding the relation between each observer elevation angle @ and the row 
numbers p in the image plane, the pixels in the image plane can be mapped into curves that use 
the same scale as Tapp _sea versus a. Finally by plotting together the previous two graphics, the 
sea Tbound can be read directly from that plot for each intersection point (Japp_sea ,a). This 
procedure can be automated by transforming the plot into an 11 by 11 matrix, for example, where 
each row corresponds to p and each column to Tbound_sea. 

To range-compensate the ship all that it is needed is to know the ship apparent black body 
temperature 7bound_ ship as a function of the apparent temperature Zapp_ship at the image plane. 
This can be done because all the ship points were assumed to have the same distance to the 
observer, so that p is constant. But a discontinuity correction is needed to adjust the ship 
temperature to the sea temperature at the sea-ship interface and complete the overall range 
compensation. 

An analysis of the range-compensated image shows that in the vertical direction the ship 


temperature has become more flattened. This is because the ship temperature is around 23°C , 
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and this temperature is a transitional temperature where 7bound-ship become higher than 
Tapp ship. For temperature less than 23 °C Tbhouna-ship is less than Tapp ship. The basis for 
that transition temperature is the magnitude of the path radiance. 

The path radiance has a great influence on temperature attenuation and it determines the 
temperature difference between Tbound and its corresponding apparent temperature in most@ 
cases. This fact is clearly seen near the horizon where the sea temperature approaches the sky 
temperature, and this will result in a range-compensated image completely different from the 
original. 

The range-compensation method proposed in this thesis has many advantages over the 
usual integral solving approach to find the original exitance of a source and range-compensate a 
pixel in the image plane. One is the reduction in the computation time because it is easier to 
interpolate than solve an integral (or to solve a matrix expression, if we change the continuous 
integral form into the discrete form). Another advantage of the proposed method is that the same 
apparent temperature matrix can be used for images taken at another time if the atmospheric 
conditions did not change dramatically, because in this case the radiance components are 
unchanged. Also the computer memory needed to solve the problem is much less in this case. 
Using a matrix form with a 5 cm’! frequency step size to solve the integral, the traditional 
approach requires the storage of vectors of 85 elements (transmittance, emissivity and radiances) 
for each row in the sea image plane. Instead, the actual method merely needs to store a matrix of 
the size of the precision desired (in this thesis an 11 by 11 matrix). But the biggest advantage is _ 
the simplicity. 

Note that the principal difference between the two methods is the direction of the 
calculation of the integral that relates the received irradiance with the sourcé exitance. The 
proposed method uses the source exitance to find the irradiance at the full range in opposition to 
the traditional method that solves for the radiant source exitance knowing the received irradiance. 
Then, both methods should give the same results if the number of rows and columns of the 
apparent temperature matrix are selected in such a way that no interpolation is needed. 

In this thesis SeaRad, Matlab, Mathcad, and other data manipulating programs were used. 
But a unique user-friendly program can be written to find the correct signal and background 
temperature and use it in predictive codes to estimate the signature at any other range. 

Although all this work was conducted for an unpolarized image, the same basic 


guidelines can also be applied for polarized cases. 
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